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Abstract 

It has been suggested about 20 years ago that the hquid close to the inner core 
boundary (ICB) is supercooled and that a sizable mushy layer has developed during 
the growth of the inner core. The morphological instability of the liquid-solid in- 
terface which usually results in the formation of a mushy zone has been intensively 
studied in metallurgy, but the freezing of the inner core occurs in very unusual con- 
ditions: the growth rate is very small, and the pressure gradient has a key role, the 
newly formed solid being hotter than the adjacent liquid. 

We investigate the linear stability of a solidification front under such conditions, 
pointing out the destabilizing role of the thermal and solutal fields, and the stabi- 
lizing role of the pressure gradient. The main consequence of the very small solidifi- 
cation rate is the importance of advective transport of solute in liquid, which tends 
to remove light solute from the vicinity of the ICB and to suppress supercooling, 
thus acting against the destabilization of the solidification front. For plausible phase 
diagrams of the core mixture, we nevertheless found that the ICB is likely to be 
morphologically unstable, and that a mushy zone might have developed at the ICB. 
The thermodynamic thickness of the resulting mushy zone can be significant, from 
~ 100 km to the entire inner core radius, depending on the phase diagram of the 
core mixture. However, such a thick mushy zone is predicted to collapse under its 
own weight, on a much smaller length scale (< 1 km). We estimate that the inter- 
dendritic spacing is probably smaller than a few tens of meter, and possibly only a 
few meters. 
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compaction, interdendritic spacing 
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1 Introduction 



As the Earth's core is g radually cool inp; down, the inner core is freezing from 
the liquid core mixture (jjacobd . Il953l ) thought to be mostly iron-nickel alloyed 
with a small quantity of lighter elements. Although th e nature and r elative 
abundances of these light elements are still controversial (jPoirierl . Il994al ). min- 
eralogical models consistent with the seismological constraints indicate that 
their concentrat ion is greater, perhaps by a f actor 4, in the liquid core than in 
the inner core (lAnderson and Ahrend . Il994j ): like most alloys, the core mix- 
ture is fractionating during the process of crystallization, the liquid core being 
slowly enriched in light elements. 



The release of light elements which results from the gradual solidification of 
the core is thought to be a ma ior source of buoyancy for driving the core 
convection and the geodynamo (IBraginskyl . Il963l : iGubbins et al.l . |200J), but 
segregation of light solute may also have dramatic co nsequences on the struc- 



ture of the solid inner core itself (IFearn et al.l . Il98ll ). It is well known from 



metallurgical experiments that solute segregation during solidification can re- 
sult in the supercooling of the liquid close to the solidification front. This 
supercooling is usually suppressed by either solidification of isolated crystals 
in a slurry layer, or by the development of a m ushy zone, where solid cells or 
dendrites coexist with a solute rich liquid (e.g. iKurz and Fishen . 119891 ). 



Loper and Roberts! (Il98ll ) and iFearn et al.l (Il98ll ) show that the conditions 
at the inner core boundary (ICB) are almost certa inly favorabl e to th e for- 
mation of either a slurry layer or a mushy layer. iFearn et al.l (Il98ll ) gave 
preference to the latter and concluded that a mushy zone of considerable 
depth, possibly extendi ng to the inner core center, must have developed at 
the ICB. More recently, IShimizu et al.l (120051 ) quantitatively studied the pos- 
sible regimes of solidification (slurry layer or mushy zone) and found that a 
dendritic regime is more probable than a slurry layer regime, because of the 
difficulty t o supply enough nucle i to fe ed a slurry layer. Some seismological 
study (e.g. ICao and Romanowica . 120041 ) also argue in favor of the existence of 
a mushy zone at the top of the inner core. 



Morsd (Il986l . l2002l ) challenged this view arguing that convective motion should 
quantitatively remove light solute from the inner core boundary, thus sup- 
pressing supercooling. From an analogy between the growth o f the inner core 
and cumulate solidification in magma chambers, iMorsd (120021 ) concluded that 
the ICB must be microscopically fiat, and rejected the idea of a mushy zone 
at the top of the inner core. Cumulates and dendritic layers have dynamics 
similar in many respects, both being reactive porous media, but their devel- 
opment occurs through two physically different processes. Cumulates form by 
sedimentation of early formed crystals from the melt, the possible porosity 
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resulting from a competition between sedimentation of crystals and solidifica- 
tion in the layer. In contrast, a dendritic layer develops as a consequence of an 
instability of an initially plane solid-liquid interface. Whether Morse's analysis 
is appropriate or not to the formation of a dendritic layer may therefore be 
questionable, but it is probably correct that advective transport of solute may 
affect the solidification regime of the ICB. 

The morphological stability of a solidification fro nt has been intensively stud- 
ied both theoret i cally and experirn entally 



Kurz and Fisherl . Il989l : iDavid . 



le.g. 



MuUins and Sekerka. 1963 



2OOII ) but, as noted by IShimizu et al.l (120051 ) 



the freezing of the inner core occurs in very unusual conditions: the growth 
rate is very small, about 6 orders of magnitude smaller than in typical labora- 
tory experiments, and the pressure gradient has a key role, the newly formed 
solid being hotter than the adjacent liquid. 

In this work, we will take into account first order effects of advective transport 
and investigate possible effects of the very slow solidification rate and of the 
pressure gradient on the stability of a solidification front (sections 2, 3 and 
4). In section 5, the thermophysical parameters and conditions relative to the 
crystallization of the inner core will be evaluated from a survey of the literature 
and from the use of thermodynamic constraints. The stability of the ICB will 
then be discussed and we shall argue that the existence of a mushy zone is 
very likely. Finally, in section 6, we will estimate various length scales relevant 
to the structure of a mushy inner core, such as its depth and interdendritic 
spacing. 



2 Supercooling in the vicinity of an interface 



When a dilute alloy is frozen from the melt, the newly formed solid usually 
differs in composition from the liquid, in a way which depends on the solu- 
bility of the solute in the solid and liquid phases. This solute redistribution is 
conveniently described by the distribution coefficient k = c^/c, where and 
c are respectively the mass fractions of solute in the solid and in the liquid. In 
the common case where k < 1, solidification results in the rejection of solute 
from the solid phase, the liquid phase being enriched in solute in the vicinity 
of the solidification front. 

As the melting temperature is lowered by the presence of solute, the liquidus 
temperature at the interface is lower than that of the liquid ahead of the 
interface, which is richer in solute. If the actual temperature gradient at the 
solidification front is smaller than the gradient of the melting temperature, 
the temperature of the liquid close to the interface is smaller than the liquidus 
temperature, i.e. the liquid is supercooled. This non-equilibrium state usually 
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does not persist and results in the formation of either a mushy zone, through 
a destabihzation of the sohdification front, or a slurry layer. 



Denoting and Tm the temperature in the liquid and the melting temperature 
respectively, the criterion for supercooling is: 



dz 



< 



dTrr 



dz 



where the subscript denotes the value at the solid-liquid interface, and where 
the z axis points toward the liquid. In the core, the melting temperature 
depends on both the solute concentration c and the pressure P. Defining the 
hquidus and Clapeyron slopes as: 



dT„ 



m, 



dc 



< and nip 



dP 



> 0, 



the supercooling criterion can be written: 



— nip Gp < rricGc 



(3) 



where On 



dT 



dz 



Gr = — 



dc 
dz 



and Gp = 

dz 



are respectively the thermal. 



chemical and pressure gradients at the solid-liquid interface. In the solidifica- 
tion conditions of t he core, the pressure gradient is negative and mp positive 
(e.g. iBoehlerl . I1993I ). thus mpGp is negative and the pressure field acts against 
supercooling. In contrast, the temperature and concentration fields both pro- 
mote supercooling, as G^ < and rricGc > 0. We can therefore predict that 
the pressure field should have a stabilizing effect on the solidification front 
whereas temperature and solutal fields should be destabilizing. 



3 The effect of advective transport on the supercooling criterion 



If motion is present in the liquid, advective transport of solute, and, to a 
smaller extent, of heat, may affect the solutal and thermal profiles near the 
solid-liquid boundary, and thus the degree of supercooling. In this section, 
we consider the effects of buoyancy-driven convection on the mean solute and 
thermal fields. Let us consider the directional solidification of a dilute binary 
alloy at constant velocity V. The solute has a concentration Coo in the bulk of 
the liquid and a chemical diffusivity in the liquid. The coordinate system is 
fixed on the moving front, the z axis pointing toward the liquid perpendicularly 
to the solidification front. The liquid has a velocity u. 

In the melt, the solute concentration c is given by the equation of conservation 
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in the bulk: 



with boundary conditions at infinity: 



at z 



oo, 



(4a) 



(4b) 



and at the interface z ■ 
the sohdification front: 



0, given by the equation of conservation of solute at 



[Co - Co J 



V 



co{k - V 



V 



(4c) 



where Cq and Cq are the concentration in the solid and liquid phases at the 
interface. If advective transport is negligible, a solute boundary layer of thick- 
ness Sc = Dc/V will build up, and after a transient stage during which the 
solute concentration in the solid at the solidification front gradually increases 
from kcoo to Coo, the system reaches a steady state where the concentration of 
solute c{z) in the liquid is 



C{z) = Coo + Co 



k 



exp {-z/6c) 



(5) 



(e.g. iKurz and Fisherl . Il989l : iDavid . l200ll ). The concentration in the liquid at 

the interface is Cq = Coo/k and the chemical gradient at the interface is 



k-lV 
' k 'K 



(6) 



However, if the solidification is very slow, as it is at the ICB, advection may 
affect strongly the concentration gradient at the interface. With a solidification 
velocity of 10~^^ m.s~^ and a chemical diffusivity of order 10~^ m^.s~^, the 



chemical diffusive boundary layer thickness is 5c 
the solutal Rayleigh number Ra** as: 



Ra^ 



DJV ~ 100 m. Defining 



(7) 



where /? ~ 1 (iGubbins et al.l . |2004| ) is the compositi onal expansion coefficient, 
q = 4 .4 m.s~^ is the magnitude of gravity at the ICB (IDziewonski and Anderson 
198ll ) and u ~ 10~^ m^s~^ is the kinematic viscosity of the liquid, we found 
Ra*^ ~ 10^^. This is much larger than the critical Rayleigh number for the 
Rayleigh- Taylor instability. The purely diffusive boundary layer is therefore 
subject to convective instabilities, which remove light solute from the vicinity 
of the interface much more efficiently than diffusion alone. 

The purely diffusive concentration profile described by equation ([5]) is there- 
fore not relevant to the case of the Earth's inner core solidification, and effects 
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of advective solute transport on the mean concentration profile must be inves- 
tigated. In the case of purely diffusive solute transport, Cq is equal to c^olk. 
This is the maximum allowed value as a greater value will imply a solute 
concentration in the solid greater than Cqo- If convective motions are present, 
advective transport results in the depletion of solute close to the solidification 
front, and the concentration in the liquid at the boundary will be closer to Coo 
as convection is more vigorous, whereas concentration in the solid will tend to 
kCoQ. Because solute conservation implies that the concentration in the solid 
cannot exceed Cqo and must be greater than /cCqo, the concentration gradient 
Gc at the interface is bounded by: 

Coo(fc - 1)-^ < Gc < Co ^ 



the upper bound being the value of the chemical gradient in the purely dif- 
fusive case. This gradient tends toward its lower bound as the efficiency of 
advective transport is enhanced; in what follows, we will consider that the 
actual chemical gradient can be approximated by its lower bound: 

Gc^Coo(A;-l)^, (9) 

which should be compared with equation If is small, the convective 
chemical gradient can be considerably smaller than the diffusive gradient. 

With a no-slip boundary condition, u{x,y,0) = for all {x,y), equation (Hal) 
taken at 2; = gives the second derivative of the mean concentration profile 
at the interface: 



Dr dz 



-yGc (10) 

Or. 





In what follows, the vertical variations of the mean concentration field in the 
vicinity of the interface will be described by a second order Taylor expansion: 

1 z"^ 

c{z) =co + G,z-G,-—. (11) 

Or I 



Because the thermal diffusivity is much larger than the solute diffusivity Z?c; 
advection affects the solute flow much more strongly than the thermal flow, 
and we will ignore the effects of advection on the thermal field, i.e. the mean 
temperature is solution of the equation of diffusion in the liquid and in the 
solid: 

- V-^ = D'V^T' (12) 



dz 

arps 

= D^V2T^ (13) 

oz 
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where and are the thermal diffusivities in the hquid and sohd respec- 
tively. The thermal gradient at the interface is given by the heat balance at 
the solid-liquid interface: 



(14) 



where is latent heat per unit volume, k'^ and are the thermal conduc- 
tivities in the solid and liquid respectively and and are the gradient of 
temperature at the interface in the solid and liquid respectively. The temper- 
ature field can be expanded as 



T(z) 



If! 



(15) 



where 5t = Dt/V. 



Shimizu et al.l (120051) did not take into account effects of convection on the 
thickness of the boundary layer, and assumed it to be given by the diffusive 
boundary layer thickness 6c = Dc/V. They study the possibility that a fraction 
of the solidification occurs in the supercooled zone, whose thickness has been 
taken to be the diffusive boundary layer thickness, that is ~ 100 m. With 
this value, they found that it is unlikely that a significant fraction of the 
solidification occurs in a slurry layer, because of the difficulty of supplying 
continuously enough nuclei. Taking into account advective transport leads 
to a much thinner supercooled zone, with smaller supercooling, where the 
amount of cr ystals solidified is pr obably much smaller than expected from 
the model of IShimizu et al.l (120051 ). We will therefore consider that all the 
solidification occurs at the solid-liquid interface, and that growth of crystals 
in the supercooled zone is negligible. 



4 Linear stability analysis of the growth interface of a binary alloy 
in a pressure gradient 



4.1 Solutal and thermal instabilities 



Considerations of section 13] allow us to determine a mean state for the solute 
and thermal fields, hereafter noted c and T, respectively described by equa- 
tions flTTl) and (IT^ : the pressure field is taken to be hydrostatic, the pressure 
gradient being Gp = —p^g. 

In this section we consider the stability of the planar interface against in- 
finitesimal perturbations of the mean fields. The interface is not planar any- 
more, but has an infinitesimal topography h{x, y, t) which temporal evolution 
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we study. The solute and thermal fields can be expressed as the sum of the 
mean field and infinitesimal disturbances : c{x, y, z, t) = c{z) + c{x, y, z, t) and 
T{x, y, z, t) = T{z) + T{x, y, z, t). The linear stability an alysis we propose here 



is sim ilar in principle to Mullins and Sekerka's analysis (IMullins and Sekerka 



19631 ). The two main differences concern the formulation of the basic state, 
which is considered here to be altered by convective motions, and the depen- 
dence of the melting temperature on pressure. 

The thermal and solutal fields must satisfy: 

(i) the equations of conservation of solute and heat in the liquid and solid 
phases: 

I - ^1 - V^c. (16a) 

f - - fl^v^r', ,16b, 

Diffusion of solute in the solid is neglected. 

(ii) the boundary conditions at the interface. At thermodynamic equilibrium, 
the temperature at the interface must obey the Gibbs-Thomson relation: 

= T^ + TH + mpGih + rUcC, (16d) 

which states that the temperature at the interface is equal to the melting tem- 
perature of the curved interface, when variations with pressure, concentration 
and interface curvature are taken into account. Tm is the solvent melting tem- 
perature of a flat interface at z = 0, ~ V^/i is the interface curvature and 
r = T^'-f/L^ is the Gibbs coefficient, where 7 is the liquid-solid interfacial 
energy. 

In addition, the thermal and solutal flelds must respectively satisfy the heat 
balance and the solute conservation at the interface: 

L.Vnix, y) = [k'VT - k'VT') ■ n(x, y), (16e) 

{c'' - c)vn{x,y) = D^Vc- n{x,y), (16f) 

where Vn{x, ~ + dh/ dt is the speed of the front and n is the unit normal 
vector pointing toward the liquid. 

(iii) boundary conditions at inflnity: the perturbations must decay to zero at 
inflnity (but see discussion below). 

Equations (fT6|) lead to a system in perturbation quantities h, c, and T'^ 
which is linearized, allowing to seek solutions for the perturbed flelds under 
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the following normal mode form: 



c 



exp{u!t + ikxX + ikyu) 



(17) 



where u is the growth rate of the disturbance and kx and ky are the wave num- 
bers along the interface in the direction x and y. From equations f|T6k .b.c) and 
the boundary condition (iii), Ci{z), T({z) and Tl{z) are found to be propor- 
tional to exp(/5z), exp{P^z) and exp{P''z) respectively, where /3, and 
are 




> 0, 



V 



V_ 



Using expression (ITTl) in the perturbation equations leads to a set of homoge- 
neous, linear equations in hi, ci, Tf and which have nontrivial solutions if 
the following dispersion equation is satisfied: 



CO ={ anicGc — a^Gx "~ "''G't + mpGp — Vk\ 

— ■ 777 ( T ^ 
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where: 



a 



a 



a 



1.2 _|_ 1.2 



P + V/Dc 
(3+{l-k)V/D, 



K 



e [0 1]. 



Stability of the interface against infinitesimal perturbations depends on the 
sign of u: the solidification front is stable if u is negative for all wave numbers 
whereas it is unstable if u is positive for any wave number. As K'^j3^ — K^l3^ > 0, 
PDc/V + 1 — A; < and rricGc > 0, the denominator is always positive and so 
the sign of u depends only on the sign of the numerator, which is a combination 
of temperature, concentration and pressure gradient weighted by wave number 
dependent functions. 

Convection has been taken into account when estimating the mean chemical 
gradient, but to be rigorous, we should have considered the effe cts of convec 



tion o n t he perturbed fields as well. We might have imposed, as ICoriell et al. 



( 119761 ) or lFavier and RouzaudI (Il983l ) did, that perturbations of the solute field 
must decay to zero on a finite length, imposed by convection, but numerical 
calculations not presented here show that this does not significantly affect the 
prediction of the stability of the interface. This can be understood as follows. 
With the additional hypothesis of a small solutal Peclet number V/ k^Dc (that 
is, for instability wavelengths small compared to 21^0 c/V ~ 300 m, an hy- 
pothesis which will be shown in section 5 to be well verified), (3 is very close 
to —kh, and a is very close to one because kh ^ V/ Dc- If taken into account, 
advection of the perturbed solutal field would decrease the decay length of 
the perturbations, and so increase j3. This would make a to be even closer to 
one, but will hardly affect the value of the chemical term in the numerator. 
Effects on the stability limit and on the critical wavelength would therefore 
be insignificant. However, an increase of (3 would result in a decrease of the 
denominator, and hence in an increase of u: the instability growth rate may be 
underestimated by equation f[T5]) . but this point will not weaken our analysis 
of the ICB solidification regime. 

The dispersion relation confirms the qualitative differences between constant 
pressure solidification experiments and the crystallization of the inner core, 
as seen from the criterion of supercooling. In contrast with usual solidifica- 
tion, temperature gradients are negative and destabilizing, as a is positive. As 
mpGp is negative, the pressure field stabilizes the interface against topography 
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perturbations. 



The differences between usual solidification and crystallization in a pressure 
field are illustrated in a In ^ versus In c diagram (figure [1]) constructed from 
equation ([T5]) . The dashed curve is the neutral curve for the purely solutal sta- 
bility problem: temperature gradients are taken to be constants and rejection 
of latent heat at the interface is neglected. In this limiting case, the neutral 



curve is very similar to that of constant pressure solidification (e.g. [Davis 



2OOII ): the curve possesses two asymptotic straight lines of slopes -1 and +1 



which correspond to the constitutional supercooling limit and to the absolute 
stability limit respectively, where the stabilizing effect of the surface tension 
becomes dominant. In constant pressure solidification, allowing perturbations 
of the thermal field slightly stabilizes the in terface for small c, and flattens 
the nose of the marginal curve (jPavisl . 1200 ll ). In the system considered here, 
the thermal gradient is destabilizing, and even when solidifying a pure melt 
(no solutal destabilization), the solidification front may become unstable if 
the solidification velocity is high enough. The neutral curve of the complete 
thermo-solutal problem is shown as a solid line in figure [H For low Coo, the 
instability is mainly thermally driven and the critical solidification velocity 
is independent of the concentration. As for the purely solutal destabilization, 
the system reaches an absolute stability limit for very high solidification rates, 
with the neutral curve being independent of Coo at small concentration. 



Effects of advection on the stability limit is illustrated in figure [H which rep- 
resent neutral curves obtained with two different expressions for the chemical 
gradient at the interface. The dashed curve has been obtained with the upper 
bound of Gc, in inequality ([8]), and therefore corresponds to a purely diffusive 
solutal transport. The solid curve has been obtained with the lower bound of 
Gci which corresponds to the case of maximum advective transport of solute. If 
the segregation coefficient k is small, convection has a considerable stabilizing 
effect: at a given solute concentration c, the critical solidification velocity can 
be as much as an order of magnitude greater than in the case of no convection. 



4- .2 Damping of instabilities by solid deformation 



The only effect of the pressure gradient which has been considered so far is the 
pressure dependence of the melting temperature. Yet, because of the differ- 
ence of density Ap between the solid and the liquid, a surface topography may 
induced horizontal pressure gradients which may tend to flatten the interface, 
therefore acting against its destabilization. To quantify the possible impor- 
tance of solid flow on the interface stability, we will estimate the timescale 
of isostatic adjustment, assuming that the viscous deformation is driven by a 
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Fig. 1. Neutral curves for linear morphological stability, constructed from equation 
([18]). The dashed line is the neutral curve for the purely solutal stability problem. 
The solid line is the neutral curve for the thermo-solutal stability problem. 

balance between the pressure gradient and the viscous force: 

= -Vp + r/V^u, (19) 

where rj is the solid state viscosity. Let us assume that the interface has a 
topography of amplitude h which varies on a length scale A. We are still 
dealing with linear stability analysis and infinitesimal perturbations, so that, 
as A is finite, h <^ X. The finite amplitude case, i.e. dendrites compaction, 
will be considered in section 6.1.2. Because h X, the deformation induced 
by the topography must be accommodated in depth, on a length-scale ~ A, 
and the horizontal and vertical velocity u and w must be of the same order 
of magnitude. The horizontal pressure gradient is of order Apgh/X, V^m is of 
order m/A^, so that: 



u 



w 



-hX. 



V 



(20) 



The timescale of isostatic adjustment can be defined as the ratio of the topog- 
raphy to the vertical velocity r = h/w, and is equal to: 



T ~ 



V 

ApgX' 



(21) 
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Fig. 2. Neutral curves for purely diffusive transport of solute (dashed line) and for 
advective transport (solid line), rric = —10^ K, A; = 0.2. 

In the limit of h <^ X which is considered here, r has a finite value even 
for an infinitesimal amplitude h of the topography. This process is therefore 
relevant when dealing with linear stability analysis. We will however consider 
that this effect is negligible compared to the destabilizing effect of the thermal 
and compositional gradient, an assumption which will be justified a posteriori 
in section [5l2l where r will be compared to the timescale l/u of the growth of 
an instability. 



5 Morphological stability of the ICB 



5. 1 Thermo-physical parameters, composition, growth rate 



As discussed in section 2, the degree of supercooling and the stability of the 
boundary depend strongly on the abundance of solute and on the phase di- 
agram of the crystallizing alloy, i.e., the liquidus slope and the segregation 
coefficient. The nat ure and relativ e abundances of light elements in the core 
are still uncertain (jPoirierl . Il994al ). but recent studies seem to show a pref- 
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erence for O, S and Si as major ligh t elements (IRingwood and Hibberson 



19911 : IStixrude et all . 119971 : lAlfe et all . 12002 : iRubie et all . l2004f ). Very little 



is known about the phase diagrams of the candidate alloys and it is not 



ior ( 


Williams and Jeanloa, 


1990 


Sherman , 


1995; 


Boehlei . 


19961). 



phase diagram of the iron-major light elements system is that the outer core 
is richer in light solute than the inner core. Major light components must 
therefore have a segregation coefficient k smaller than one and this necessarily 
imposes that their liquidus slopes rric are negative near the Fe-rich end. 



The liquidus slope and solute concentration of interest here are that of the 
fractionating elements. Elements which do not frac tionate du r ing th e process 



of solidification, such as Ni and perhaps S and Si (lAlfe et al.l . |2002| ). will not 
create chemical heterogeneity and hence will not contribute to supercooling 
and radial variations in melting temperature. 



Assuming ideal mixing, a crude estimate of the hquidus slope of the core 
mixture at ICB pressure and temperature may be provided by the van't Hoff 
relation (j Chalmers! . Il964j ). We consider here the effect of alloying Fe with a 
single light element, of mole fraction x, on the melting temperature of the 
mixture. Chemical equilibrium between two multicomponents phases requires 
equality of the chemical potentials of each component in the two phases. In 
particular, the chemical potentials of the solvent, here Fe, must be equal in 
the liquid and in the solid: /Xpg = fip^. In an ideal solid or liquid solution, the 
chemical potential of a component i is expressed as fii = fi° + RTliaxi, where 
Xi is the mole fraction of component i and /i° is the chemical potential of pure 
i. Equality of the chemical potentials of iron in the liquid and solid phases 
then requires that 



fi'^l + RTlnx' 



.s 
Fe 



+ i?T In 4„ 



(22) 



which can be rewritten, using the fact that the mole fraction of the solute is 
X = 1 — xpe and that x'^ = kx^, as 



,o s _ , ,o £ 2 



T 



Rln 



1 — kx^ 



(23) 



Taking the derivative of equation (1231) with respect to T, and using the Gibbs- 
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Helmoltz relation then gives 



5 A 1 — ' 




(24) 



"'Fc "'Fe 



where /ipf and h^^ are the molar enthalpy of solid and liquid Fe respectively, 
L = {hpl — hp^)/Mpe is the latent heat of pure iron and Mpe is the atomic 
weight of iron. Assuming k to be constant, for a dilute solution {x^ <^ 1), we 
obtain 





which together with 0241) yields, as mc can be written as dTm/dx^ 



(25) 



p'T'o 2 



(6±3)(A;-1) X lO^K 



(26) 



where is the melting temperature of pure iron at ICB pressure and ttIc 
is given in Kelvin per atomic fraction; values and incertitudes of and L 
are from table [H This relation means that for an ideal solution, the liquidus 
slope and the segregation coefficient are related by parameters which are inde- 
pendent of the nature of the alloying element. For a non-partitioning element 
(fc = 1), the ideal liquidus slope is equal to zero, whereas a highly partitioning 
element {k small compared to 1) will have a high liquidus slope. Through this 
relation, constraints on k may provides informations on the liquidus slope of 
the core mixture. 



From estimates of the volur ne change during melti n g and of the density jump 
at the inner core boundary, [Anderson and Ahrend (119941 ) estimated the ratio 
of light elements in the outer core to that in the inner core to be approximately 
4 to 1. This gives a global segregation factor k = 0.25 and, by equation fl26|) . 
a liquidus slope rric — —4.5 ± 2.5 x 10^ K. Note that the segregation factor 
estimated here is an effective segregation factor, which is higher than the 
thermodynamic one, and should give a lower bound of \mc\- This estimate may 
be appropriate if there is only one light element in the core, but finer estimates 
are needed if there are several light elements of comparable c oncentrations. 
As an example, the ab initio simulations of lAlfe et al.l (120021 ) suggest that 
the outer core may be composed of ~ 10 mole% of S a nd /or Si and ~ 8 
mole% (~ 2 wt.%) of O. According to lAlfe et all (j2002f ). S and Si do not 
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significantly fractionate and tlie density jump at tlie ICB may be accounted 
for by fractionation of oxygen alone, whose segregation coefficient has been 
estimated to be 0.02. With this value, rric tends to its = bound rric — 
—6 ± 3 X 10^ K. In what follows, two chemical models of the core will be 
considered: one with a single light element of concentration ^ 10 wt.%, and 
the other with only one fractionating light element (but several light elements), 
oxygen, of concentration ~ 2 wt.%. 



Our r elati on is similar to the one d erived by lAlfe et al.l (120021 ). [Stevenson 
( 1l98ll ) and [Anderson and Dubai (119971 ) estimated the melting point depression 
by assuming equilibrium between a pure solid and an alloyed liquid, and there- 
fore obtained the upper bound of our estimate (i.e. our k = value). Those 
theoretical estimates are in poor agreement with experi ments. Experimental 
resul t s on the melti ng temperature of the Fe-0 system (IKnittle and Jeanlozl . 
199ll : iBoehlerl . Il993l ) predict a small melting temperature depression, and per- 
haps a^o]id2Sohjtion__bel^ case of the Fe-S system if more controver- 
sial. Williams and Jeanloz (Il990f ) results suggest that the eutectic behavior 
persists at high pressure, and predict a significant melting point depression 
(rrir, is of order —5 x 10^ K at core-mantle boundary pressure). In contrast 



Boehlerl (119961 ) found that the Fe-FeS eutectic melting depression becomes 
much smaller at high pressures, and conclude that this supports the possibil- 
ity of solid-solution between Fe and FeS at core pressures. To our knowledge, 
there is no experimental work at this pressure range dealing with other can- 
didate alloys. In the present work, values of tjIc between —10^ K and —10"^ K 
will be considered. 



The interfacial energy of iron at ICB conditions can be deduced from estimates 
of the latent heat of crystallization, because those two parameters both derive 
from the difference of atoms bonds energy between the solid and liquid phase. 
The interfacial energy per atom 7^ can be c alculated to be 1/4 of the atomic 
latent heat La for a fiat close-packed surface (j Chalmers! . Il964j ) . Estimates of the 
l atent h e at of iron at ICB conditions range from 600kJ kg~^ to 1200kJkg~ ^ 



( Poirien. Il994bl : lAnderson and Dubai . 119971 : iLaio et al.l . I2OOOI : IVocadlo et al. 



2003al ) and from these values, we estimate the interfacial energy per unit area 
7 to be 0.4 ± 0.2 J m~^, which can b e compared to the 0.204 J.m~^ value at 
standard conditions (jChalmerd . Il964l ). 



Buffett et al.l (119921 ) proposed an analytical model of growth of the inner core 
and found that the radius increase at leading order as the root square of time, 
r = He \Jt/a, where r^^ is the present radius of the inner core and a is its age. 
The present solidific a tion velocity is then V = r^ J'^a. If the inner core is young 
{e.g. iLabrosse et al.l . I2OOII : iNimmo et al.l . 120041 ). i.e. a ~ 1 Ga, V is found to 



be 2 X 1 m.s ^. On the other hand, if the inner core nucleated around 3 Ga 
ago, a s IChristensen and Tilgnerl (120041 ) claim, V could be of order 6 x 10~^^ 



m.s 



-1 



WenI (I2OO6I ) observed a temporal change of travel time of the PKiKP 
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phase between the two events of an earthquake doublet, indicating a locahzed 
change of the inner core radius of about 1 km in ten years. This observation 
may be interpreted as reflecting episodic growth of the inne r core, coupled 
with non stationary convection in the outer core (jWeru . |2006| ). The resulting 
instantaneous solidification velocity is ~ 10 km/10 years ~ 10~^ m.s~^, which 
is much higher than the mean solidification velocity estimated from models of 
the core thermal history. 



Other parameters used in this study, with values currently found in the liter- 
ature, are listed in table [H 



Table 1 

Parameters used in this study. 



Clapeyron slope 


mp 


~ lO-S KPa-i « 


Liquidus slope 




-10^ to -10^ K ^ 


Light elements concentration 


Coo 


2 to 10 wt.% 


Growth rate of the inner core 


V 


6 X 10-^2 to 2 X 10 


Temperature at the ICB 




5000 to 6000 K 


Specific heat at constant pressure 


CP 


860 Jkg-iR-i 


Latent heat of crystallization 


L 


600 to 1200 kJkg 


Thermal conductivity in the liquid 




63 Wm-iR-i <^ 


Thermal conductivity in the solid 




79 Wm-^K-i 


Thermal diffusivity in the liquid 




6 X 10-^ m^s-^ 


Thermal diffusivity in the solid 




7 X 10-6 m^s-i 


Chemical diffusivity 


Dc 


10-9 m^s-^ " 


Viscosity of the inner core 


V 


10^*^ to 10^1 Pas 


Liquid-solid interfacial energy 


7 


0.4 ±0.2 Jm-2 ^ 


Gibbs-Thomson coefficient 


r 


~ 2 X 10-^ K.m ^ 



-1 a 



" from I Anderson and Dubai (|1997| ) 
^ see text, 
from 



from 



Poirieil (Il994bl'). 



Stacev and Anderson ( 200ll ) 
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5.2 Supercooling and stability analysis 



The liquid at the ICB is supercooled if condition ([3]) is satisfied. This criterion 
can be rewritten, using the heat balance at the interface (equation (1141) ). as: 

L V 

-^ + m^Gc-—G?Y>pgm?. (27) 

Depending on the age of the inner core, on the value of the thermal conduc- 
tivity and on the hypothetical presence of radioactive elem ents, the therma l 



gradient in the inner core at the ICB may vary widely (lYukutakd . Il998l ). 
However, the thermal term L^V^//t^, which is the contribution to supercooling 
from the heat released by crystallization, is high enough to balance the pres- 
sure term if the solidification velocity is higher than ~ 5 x 10~^^ m.s~^. This 
is roughly equal to the lowest estimates of V , and, as the term — is posi- 
tive, it is thus likely that thermal terms alone are sufficiently high to ensure 
supercooling. The chemical term rUcGc is probably much higher, at least one 
order of magnitude greater than the pressure term if rUc is as small as —10^ 
K and three order of magnitude greater if nic = —10^ K, and the conclusion 
that the vicinity of the ICB is indeed supercooled seems to be inescapable. 

We now use the dispersion equation (fT8|) to investigate more quantitatively 
the stability of an initially plane solid-liquid interface in the actual conditions 
of the ICB. As discussed in section 3, we choose to take the lower bound of the 
chemical gradient, given by equation (Q, so that the effect of convection on 
the basic state may be overestimated rather than underestimated. We tested 
several values of each parameter within their uncertainty ranges and found 
that for all plausible sets of values, the ICB is unstable. The liquidus slope is by 
far the most critical parameter. Because the chemical gradient is proportional 
to nic, the two orders of magnitude uncertainty on nic propagate directly 
in the uncertainties on the location of the marginal stability curve. Neutral 
curves for linear stability are plotted in figure [3l with itLc = —10^ K, —10^ 
K and —10^ K, = K.m~^ and k = 0.25. Despite the strong dependence 
on of the neutral curve location, the solidification velocity is two orders 
of magnitude greater than the critical velocity for a hquidus slope as small 
as —10^ K, and three or four orders of magnitude greater than the critical 
velocity for a liquidus slope of —10^ K or —10^ K. If the growth of the inner 
core is episodic, with a much higher instantaneous solidification velocity, the 
ICB would be even more unstable. 



In figure HI the growth rate of infinitesimal perturbations has been plotted 
against wave length, at the conditions of the ICB, for three different values 
of nic. The corresponding time-scales range from a few years to about 300 
years. We may now estimate the timescale of isostatic adjustment, given by 
equation (ETi) . and compare it to the timescale of instability growth. The 
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Fig. 3. Neutral curves for linear morphological stability, with mc = —10^ K (dotted 
line), rric = —10^ K (dashed line) and rric = —10^ K (solid line). G?^ = K.m~^ 
and k = 0.25. The straight lines cross is the location of the ICB in the stability 
diagram. Uncertainties in V result from uncertainties in the age of the inner core. 
The value of Cqo depends on the chemical model chosen: Cqo — 1 wt.% if there i s 
only one dominant light element, and Cqo — 2 wt.% if the model of lAlfe et al.l (l2002l ) 
is adopted. 



solid inner core viscosity is p oorly constrai n ed: estimates r ange from 10^^ 
Pa.s to more than 10^^ Pa.s fiBuffettl . Il997l : lYoshida et all . 1 19961 ). A lower 
estimate of the viscous flow timescale may be given by taken r] = 10^^ Pa.s. 



With Ap = 600 kg.m ^, g = 4.4 m.s ^ (IDziewonski and Andersonl . Il98ll ). the 



isostatic adjustment timescale is of order 10 s ~ 3 x 10 years for A ~ 1 cm, 
about four or six orders of magnitude greater than the timescale of instability 
growth: isostatic adjustment will not delay the insta bility. Our results are 
consistent with the conclusion of IShimizu et al.l (120051 ) that the timescale of 
dendrites growth is very short compared to the timescale of inner core growth, 
suggesting that a mushy layer will indeed form at the ICB. 

Results have been presented here for parameters (growth rate, thermal and 
pressure gradient) inferred for the present state of the inner core, but when 
dealing with the internal structure of the inner core, it is of equal interest 
to investigate what was the solidification regime during its past history. As 
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X 10 




Fig. 4. Growth rate of infinitesimal perturbations against wave length, at the con- 
ditions of the ICB, for rUc = - W^ K (solid line), -10^ K (dashed line) and -10^ K 
(dash-dotted line), k = 0.02 and = 600 kJ.kg-^ 



explained before, the solidification regime depends primarily upon the ther- 
mal, solutal and pressure gradients at the interface, the thermal and solutal 
fields being destabilizing whereas the pressure field is stabilizing. Because the 
solidification velocity was most certainly greater in the past than it is today, 
the rates of release of heat and solute were also greater; thermal and solutal 
gradients were therefore steeper (more destabilizing). In addition, the liquidus 
slope mpGp is less steep at deeper depth (because the gravity field, and hence 
the pressure gradients decrease to zero at the center of the Earth), and there- 
fore less stabilizing. At first order, all terms seem to act in the same way, and 
it is then likely that the solidification of the inner core has been dendritic for 
most of its history. 
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6 Length scales of the mush 



6.1 Vertical length scales 



6.1.1 Thermodynamic depth of the mushy zone 

In the laboratory, and in metallurgical applications, the mushy zone depth is 
typically a few centimeters. However, because of the very small temperature 
gradient in the inner core, and of the additional effect of pressure, the tem- 
perature in the inner core remains close to the melting temperature, which 
suggests that liquid enr iched in solute may remain thermodynamically stable 



at considerable depths (IFearn et al.l . Il98ll ) . 



While freezing occurs in the mushy zone, interdendritic melt is further enriched 
in solute, lowering its melting temperature further. An enriched liquid phase 
can coexist with the solid phase as long as the actual temperature is above 
the melting temperature. The depth ^xh at which a liquid of concentration 
c = Co + Ac can be in thermodynamic equilibrium with the surrounding solid 
phase is then given by equating the actual and melting temperatures: 

T(M=^m(M- (28) 

To a good approximation, the acceleration of gravity is linear in r in the inner 
core, so that pressure is quadratic in radius. Assuming for convenience the 
temperature profile to be quadratic as well (which is a not so bad approxima- 
tion if the cooling rate is approximately constant within the inner core), the 
following expression for ^Th can be found: 



where AB ~ 150 to 200 K is the difference between the actual temperature and 
the r nelting ternperatu res (at outer core composition) at the center of the inner 



core (jYukutakd . Il998l ) . The precise form of the temperature profile (quadratic 



or not) is not of great importance for our order of magnitude estimates. 

The mushy zone can extend to the center of the inner core (5Th/''"ic = 1) if the 
variation of melting temperature with pressure is compensated by |mc| Ac. The 
maximum allowable light elements concentration is the concentration at the 
eutectic (if it exists), and the maximum value ofJrn^JAc is the refore |mc|(cE — 



Co). Taking ce = 25 at.% as a plausible value (jStevensoru . Il98ll ) . we found from 
equation (!29l) that (5Th = ^ic if I'^cl is greater than 1.3 x 10^ K. On the other 
hand, if |mc| is as small as 10^ K, the resulting depth of liquid thermodynamic 
equilibrium is of the order of 300 km. 
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6.1.2 A collapsing mushy zone? 



As noted by lLopeii (jl983l ). a fundamental observation is that the ICB appears 
to be sharp on seismic wavelength scale (~ 10 km), which means that the 
solid fraction at the top of the inner core must become significant within a 
few kilometers in depth, a fact that seems at first view hard to reconcile with 



the presence of a hundred kilometers thick mushy zone. iLoperl ( 119831 ) argued 
that convective motions in the mush may cause the dendrites to thicken and 
calculated that the solid fraction may become of order one at ~ 300 m below 
the top of the mushy zone, thus explaining the seismic observations. 

In addition to this process, the solid mass fraction may also increase in depth 
because of gradual collapse of the dendrites under their own weight. The phys- 
ical process at stake is quite similar to the isostatic adjustment considered in 
section HI the density difference between the solid and the liquid implies hori- 
zontal pressure gradients which increase with depth in the mushy zone, making 
the dendrites to broaden at their base by viscous solid flow. The timescale of 
dendrites widening may be estimated by considering equilibrium between the 
horizontal pressure gradient and the viscous force in the solid (equation (IT^ ). 
as was done is section |H Here the horizontal length scale A, say the diameter 
of a dendrite, is supposed to be small compared to the height h of a. dendrite. 
As before, the horizontal velocity at depth h is of order: 

u ~ ^hX, (30) 
V 

but the timescale of interest is now the ratio of horizontal length scale to the 
horizontal velocity Tc = X/u, so that the timescale of dendrite widening - or 
compaction - is: 

- " (31) 



Apgh' 

which is inversely proportional to the height of the dendrite. As long as the 
timescale of dendrites growth is small compared to r^, the depth of the mushy 
zone increases. When increasing the height of the dendrites, the timescale of 
dendrites widening decreases, and a stationary state is eventually reached, 
where the growth of the dendrites by solidification at their tips is balanced 
by the collapse of the mushy zone. This stationary depth is then found by 
equating to the timescale of dendrites growth, which is h/V when the system 
is stationary. This yields: 

" ^ /£■ <^^' 

Uncertainties in h come mostly from uncertainties in the viscosity. With rj = 
10^^ Pa.s, h is found to b e a few meters , a ve ry small value, whereas t] = 10^^ 



Pa.s leads to ~ 1 km (ISumita et al.1 . Il996l ). Although estimates of h span 



more than two orders of magnitude, they are in any case significantly smaller 
than the thermodynamic estimate of the depth of the mushy zone. This means 
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that as a result of compaction, the sohd fraction may become of 0(1) at a few h 
in depth, well before the thermodynamic limit of liquid equilibrium is reached. 



6.2 Horizontal length scale: interdendritic spacing 



Interdendritic spacing Ai (the distance between two dendrites tips) is a length 
scale of considerable interest for the structure of the inner core. Because 



each columnar crystal is usually made of many dendrites (IKurz and Fisherl . 
19891 ). the interdendritic spacing gives a lower bound of the grain size, and 
hence may gives hints on the deformation mechanism and viscosity relevant 
to the possible viscous deformation of the inner core. In addition, the pri- 
mary dendrite spacing is needed for e stimating the permeability of the mush, 
which is roughly proportional to (IBergman and Fearnl . Il994l ). Permeabil- 
ity is a necessary paramet e r for the study of the mush hydrodynamics (e.g. 
Bergman and Fearnl . Il994l; iLe Bars and Worsten . |2006| ). or for quantitative 



compaction models (Sumi ta et al. . 19961 ). 



We make the assumption that dendrites have an axisymmetric shape described 
by a function g{z), which is linked to the solid fraction f{z) by the relation: 



A(Ai) 



where A{Xi) is the horizontal surface area occupied by a dendrite. In a cubic 
dendritic array, A{Xi) is simply A^, whereas in a hexagonal dendritic array, 
A{Xi) = 3/2tan(7r/3)A2 - O.SGA^. 



/ and g^ are proportional through a constant involving Ai, so that if / is 
known, an additional constraint on g is enough to determine Ai. This con- 
straint can be give n by Langer and Miiller-Krumbhaar's theory of dendrite 
tip radius selection (ILanger and Miiller-Krumbhaaii . 119771 ). in which the den- 
drite tip radius R is equal to the shortest wavelength A, for which the interface 
is unstable. This theory is in very good agreement wi th experiments, and we 
will follow it in the present work. It can be shown (see iKurz and Fisherl . Il989l ) 
that at low solutal Peclet number (Pec = VR/D ^ 1), a condition well sat- 
isfied here, the chemical gradient Gc at the dendrite tip is equal to the lower 
bound of the chemical gradient we used in our stability analysis (equation IQ). 
Estimates of Aj from our stability analysis thus directly give the adequate R. 
A good approximation for R at low solutal and thermal Peclets numbers can 
be derived from equation (ITSl) . noting that the chemical ter m in the numerator 
of (flSll is much bigger than the thermal and pressure terms (IKurz and Fisherl . 
1989p 



R ~ 27r. 



YD 



mcCoo{k - 1)V 



(34) 
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R ranges from ~ 1 mm if rric = —10^ K to ~ 1 cm if rric is as low as —10^ K 
(see figure H]). Uncertainties from other parameters than rric are much smaller. 



On the other hand, the dendrite tip radius R is equal to the radius of curvature 
of g taken in z = which is by definition equal to: 



R 



9"iO) 



(35) 



where R is defined to be positive. Inserting equation fl33l) into fl35|) . with the 
boundary condition /(O) = 0, yields 



or, for an hexagonal dendritic array: 



2ttR 


dl 




dz 


ICB 



(36) 



Ai 



2TrR 


df 

0.86 




oz 


ICB 



(37) 



Although simple, this relation is very general, the only assumptions made here 
being that of stationary state and that of an axisymmetric dendrite tip. If the 
appropriate assumptions on solute transport and selection of dendrites tip 
radiu s are m ade, classical models of primary dendrite spacing (ITrivedil . Il984j : 
Huntl . 1 19791 ) follow from equation (1571) . 



These models have been derived for non-convecting mushy zones where the 
solid fraction increase smoothly from zero to one at the thermodynamic depth 
of the mush, that is, df/dz\icB scales as l/(5ThQ However, seismology tells 
us that the ICB is sharp on a scale of 10 km, which means that the scale 
length of solid fraction increase at the top of the inner core is obviously not 
the thermodynamic mushy zone depth, and must be at maximum ~ 10 km. 
Therefore a minimum order of magnitude of df /dz\icB is 1/10 km = 10~^ 
m~^. This may give upper bounds for Ai; we found Ai ~ 10 m if rric = —10^ 
K and Ai ~ 30 m if rric = —10^ K. Tighter constraints may be found if we 
could estimate more precisely df /dz\icB- Compaction is not expected to affect 
significantly the dendrites spacing because it is not effective in the very top 
of the mush. The solid fraction profile in the upper few meters, from which 
the primary spacing results, is more probably controlled by convection. In his 



^ Note that, in laboratory experiments, with no pressure gradient and with a 
constant temperature gradient S^h ~ mdce — Coo)/G't, so that, as i? oc V~^^'^, 
Ai oc Grj^^'^V~^^'^. Thi s scaling law has been shown to be in very good agreement 



with experiments (e.g. Kurz and Fisher . 19891 ). 
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convective mushy zone model. iLoperl (Il983l ) found that df ldz\ic^ ~ 3 x 10~^ 
m~^. With this value and our estimates of i?, we found Ai to be 1 to 5 meters. 
In Loper's estimate, df /dz\iQ-Q is proportional to the typical velocity W of 
descending liquid close to the ICB, taken to be around 10~^ m.s~^. This value 
is quite uncertain, but as Ai is inversely proportional to the square root of 
9//9z|icB, a change of two orders of magnitude in W is needed in order to 
change the order of magnitude of Ai. As a consequence, and considering all 
sources of uncertainties, it appears difficult to have an interdendritic spacing 
of more than a few tens of meters, our preferred estimate being a few meters. If 
the growth of the inner core is episodic, the interdendritic spacing is expected 
to be smaller. As Ai oc V^^^^^ taking V = 10~^ m.s~^ instead of = 10~^^ 
m.s~^ would result in an interdendritic spacing about twenty times smaller 
than estimated above, that is, a few tens of centimeters. 



Our estimates are significantly sma ller than estimate of a few hundred meters 
from scaling laws (IBergmanl . Il998l ) . This is not surprising because those scal- 
ing laws have been derived for non-convecting (and non-compacting) mushy 
zones, where / is linearly increasing in the whole thermodynamic mushy zone. 
Assuming for a demonstrative purpose that / increases linearly from to 1 
within the thermodynamic mushy zone, whose depth is taken equal to the in- 
ner core radius, equation (1571) gives Ai ~ 100 — 300 m, in very good agreement 
with Bergman's estimate. 



Our assumption of an axisymmetric dendritic s hape is justified if the iron 



phase at inner core conditions is fee or bcc iron (jVocadlo et al.l . |2003bj), but 



might not hold if it is e-iron, as hep mate rials (e.g. ice or zinc) usually have 
plateshaped dendrites, known as platelets (IBergman et al.l . l2003l ). If iron den- 
drites are indeed platelets, estimating the inter-platelets spacing would require 
some modifications of our analysis to take into account the specific geometry 
of platelets. Although this is expected to give quantitatively different results, 
this should not alter the qualitative conclusion that convection in the mush 
reduces the interdendritic spacing, and that interdendritic spacing at ICB 
should be much smaller than suggested by classic scaling laws. 



7 Conclusion 



The morphological stability of an initially plane solidification front at the ICB 
conditions has been investigated. Despite the stabilizing effects of convection 
and of the pressure gradient, a continuous solidification implies non-zero solu- 
tal and thermal gradients at the ICB which, for plausible parameters values, 
are high enough for the interface to be destabilized. Because the conditions in 
the past were even more destabilizing, it is probable that the ICB has been 
dendritic through most of the inner core history. 
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Thermo dynamic con s iderat ions predict a very thick mushy zone which, as 



noted by lFearn et aLl (jl98ll ). could possibly extend to the center of the Earth. 
However, considerable uncertainties on the phase diagram do not allow a pre- 
cise estimate of the thermodynamic depth, which may be only a few tens of 
kilometers if the liquidus slope is small. The most superficial part of the inner 
core may be understood as a collapsing mushy zone, where both convection 
and compaction act to rapidly increase the solid fraction within a length scale 
probably smaller than 1 km, making the ICB to appear seismically sharp. 
The length scale of the thermodynamic depth of the mushy zone and the 
compaction length scale, although clearly different, are unconstrained by the 
current knowledge of the phase diagram of the core mixture and of the solid 
iron viscosity at inner core conditions. Progress in the determination of these 
parameters may greatly help the understanding of the inner core structure. 

We tried to constrain the primary dendrite spacing Ai of the mushy zone. Ai 
appears to depend on the vertical derivative of the solid fraction at the top 
of the mush which, in turn, depends on the vigor of convection in the mush. 
Here again, it is difficult to make precise and reliable estimates, but we found 
that the interdendritic spacing is most probably smaller than a few tens of 
meters, and possibly only a few meters. 

Whether or not a significant amount of melt may subsists at larg e depth re- 

main s an open question. If the mush permeability is of order A^/100 (IBergman and Fearn 



I994J ). our e stimate of Ai sugges ts permeability values higher than 10 ^ m^ 



The work of ISumita et al.l (119961 ) suggests that with such a high permeability, 
the liquid will be very efficiently removed from the inner core by compaction, 
and that the residual liquid fraction will be essentially zero. A high permeabil- 
ity in the mush does not however rule out the possibility that unconnected, 
trapped liquid pockets persist in depth. The answer to this question may de- 
pend in part on surface tension driven processes. Sintering, i.e. migration of 
grain boundaries driven by surface tension, may play an important role in the 
redistribution of the liquid phase and have consequences on the efficiency of 
compaction, in a way which will depend on the wetting properties of the melt. 
If the (unknown) dihedral angle o f light elements ric h liquid iron in contact 



with solid iron is greater than 60° ( iBulau et al.l . 119791 ) . the liquid phase is ex- 



pected to become unconnected at a given liquid fraction, thus living a residual 
liquid phase in the inner core. 
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